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Abstract. We analyze splitting algorithms for a class of two-dimensional fluid 
equations, which includes the incompressible Navier— Stokes equations and the 
surface quasi-geostrophic equation. Our main result is that the Godunov and 
Strang splitting methods converge with the expected rates provided the initial 
data are sufficiently regular. 



1. Introduction 

Let r > be a finite final time. We are interested in solutions 9: [0, T] x M 
to the generalized active scalar equation 

6*4 + M • V6' + A"6l = 0, in(0,T)xM^ a€(0,2], (1.1) 

where V is the gradient operator and A = (— A)^/^ is the fractional Laplacian 
defined through Riesz operators (see Section [2]). The divergence free velocity u is 
determined directly from through the nonlocal relation 

u = cm:\K-^e, /3e[l,2], 

where curl — V^!: denotes the spatial curl operator defined by curl(0) = {—(f>y, 4>x)- 



The general active scalar equation (1.1) seems to have appeared first in the 



mathematical literature in [3 . The general formulation encompasses a whole class 
of two-dimensional fluid equations, interpolating between the Euler/Navier-Stokes 
equations and the surface quasi-geostrophic equation. Different choices of a and 
/3 lead to different fluid equations. The most interesting (and studied) examples 
are the Navier-Stokes equations (a — /3 — 2) and the surface quasi-geostrophic 
equation (/3 = 1, < a < 2) [F; in the latter case, u = {—TZ26,TZiO) with TZi,TZ2 
denoting the usual Riesz transforms in M^. 

The quasi-geostrophic equation has recently received considerable attention from 
the mathematical community. In particular, since the two-dimensional Navier- 
Stokes equations admit smooth solutions, it has been a question if the quasi- 
geostrophic equation exhibits similar behavior. In this respect, it is common to 
distinguish between three cases of a for the geostrophic equation. When a G (1, 2), 
the dissipation term A" provides enough regularization to guarantee the existence 
of smooth solutions [B]. When a £ (0, 1), the global properties of solutions are still 
open. The remaining case (a = 1) is known as the critical case, and the global 
behavior of solutions was settled only recently (cf. [H [Ml |8] and the references 
therein) . 



In terms of the physical applicability of ( 1.1 1, it seems like the most relevant mod 



els are the critical geostrophic equation and of course the Navier-Stokes equations. 
In particular, the critical geostrophic equation has been proposed as a simplified 
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model for strongly rotating atmospheric flow. We refer the reader to [H] for more 
on the physical aspects of the model. 

We now turn to the main topic of the present paper, namely operator splitting 



algorithms for computing approximate solutions to (1.1). Generally speaking, the 
label "operator splitting" alludes to the well-known idea of constructing numerical 
methods for an intricate partial differential equation by reducing the original equa- 
tion to a succession of simpler equations, each of which can be handled by some 
efficient and tailor-made numerical method. The operator splitting approach has 
been comprehensively described in a large number of articles and books. We do 
not survey the literature here, referring the reader to the bibliography in P]. 



Regarding the Navier-Stokes equations, which is a special case of ( 1.1 ), operator 
splitting (viscous splitting) has been analyzed and applied in a great number of 
works, see, e.g., the book by Majda and Bertozzi [17]. Indeed, viscous splitting has 
been frequently utilized as a design principle for numerical methods for the Navier- 
Stokes equations, including vortex or particle methods and transport-diffusion or 
characteristic-Galerkin methods. Error estimates for Godunov and Strang viscous 
splitting algorithms have been established, e.g., in [T7], utilizing arguments that are 
different from ours and which rely on the well-known fact that there exist unique, 
global smooth solutions to the two dimensional Euler and Navier-Stokes equations. 



In this paper we apply operator splitting to separate the effects in (1.1) of the 
transport term u = curlA~^6'- V6' and the fractional diffusion term A"6'. This type 
of splitting is reasonable as it allows for specialized hyperbolic methods to be applied 
in the transport step and specialized "Fourier space" methods in the fractional 
diffusion step. The interested reader can consult [9] for further information on 
operator splitting. 

Our main contribution is that we contribute rigorous proofs of the expected 
convergence rates for operator splitting applied to the general active scalar equation 



(1.1). Our approach is inspired by the recent paper [I^ (see also [H]) on splitting 



algorithms for the KdV equation, and for that reason our results apply under the 



standing assumption that there exists a smooth solution to (1.1 ). This assumption 
is verified for the Navier-Stokes equations and the quasi-geostrophic equation with 
a > 1. It is also reasonable to expect the existence of unique smooth solution in 
the regime a g [1,2] and /3 e [1,2], cf. d [12 [HI [E] for results in that direction. 
Let us now discuss our splitting methods in more details. For this purpose, we 



write ( 1.1 1 in the form: 

et = c{9), c{9) ^ A{e) + B{e), 
B{6) = -cnriK-i^e -ve, A{e) = -K°'e. 

We will need the solution operators ^A{t, Gq) and ^B{t, do), defined as the solutions 
to the abstract differential equations: 

5t$A(i, ^o) = ^($A(t, do)), $a(0, Bo) = ^o, 
dt<fBit,eo)^B{^B{t,eo)), $a(0,(?o) = ^0- 

For local- in-time existence results for the inviscid equation 9t + u ■ V9 — (i.e., 
existence of the $b operator) when the initial data belong to Sobolev or Triebel- 
Lizorkin spaces and P = I, see [S1[71I2]; although for /3 ^ 1 such results cannot be 
found in the literature, they can be proved by properly adapting the arguments in 
[SI [3 [5] . Regarding the $a operator, the fractional diffusion equation 

z;t + A"0 = O, v{{))=Vo, 

has a solution given by v{t) = Gait)*vo, where Ga{t,x) is the fundamental solution 
in which can be expressed in terms of the Fourier transform Qa{t,C) = e^*l^l . 
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If Ga{x) denotes the inverse Fourier transform of e 1^' , then 



We refer to [TBJ [20] for further information. 

The first operator sphtting method we wiU study is known in the hterature as 
Godunov splitting. The method is defined as follows: Set 0'^ ~ 0q and sequentially 
determine approximate solutions 0", n = 1, . . . , M, satisfying 

r = $^(Ai,$B(At,r-i)) = $^(Ai) o $s(AO(0"-'). 

Formally, one can show that |j6'(nAi) — ^"ll = ©(At) in an appropriate spatial norm 
in the limit At — ^ and nAt = t. In Section [2^ we will rigorously prove this linear 
convergence rate. Specifically, we show that (for small At) 

\\0{nAt) - B^H"-- < C{T)At\\04H>'. 

for all fc > 5. The second method we will consider is known as Strang splitting: Let 
0^ = 00 and determine sequentially 

nn+l_^ /At ^ ,At 



$B(^,<I'A(At,$B( — ,r)) 
*B(^)°*A(At)o$B(^) 



1- 



Formally, one can show that the method is second order in At. That is, \\0{nAt) — 
6*" II = 0{At^) in an appropriate spatial norm in the limit At — > and nAt = t. In 
Section 12.21 we show that 



\\0{nAt) - 0^H>.-3. <Gi\\0o\\H-)At^ 
for any sufficiently small time step At > and for all k > max{5, 3a}. 



To build fully discrete numerical methods for the fluid equation ( 1.1), we have to 
replace the exact solutions operators and $b by appropriate numerical methods. 
However, we will not discuss that here. 

The paper is organized as follows: Section [2] is of an introductory nature and 
collect some results to be used later on. The convergence rate result for the Godunov 
splitting is proved in Section [3j while the Strang splitting is analyzed in Section [4j 

2. Preliminary material 
2.1. Existence and regularity results. Presently there is no complete existence 



theory for the equation (1.1). In this paper, we will assume the existence of a 
unique solution with the same regularity as the initial data. In the literature, one 
can find results confirming this assumption for some specific cases of a and l3. 



For a = (3 = 2, the equation (1.1) is the incompressible Navier-Stokes equations. 



The following result is by now classical (cf. [T5]). 

Theorem 2.1 (Navier-Stokes). Let a = /? = 2 and T > 0. If 0o e H'' , there exists 



a unique solution £ G{0,T; H ) of (1.1) with initial data 0\t=o 



The next theorem gives the well-posedness of the quasi-geostrophic equation 
when a G [1,2). The sub-critical case (a > 1) was established in |6j. Well-posedness 
in the critical case a = 1 can be found in [1] HU |H] . 

Theorem 2.2 (quasi-geostrophic). Let l3 — 1, a G [1,2), and let T > be a final 
time. Assume that 0q G for some fc > 1. Then, there exists a unique solution 



e G{0,T;H'') of (1.1) with initial data 0\t=a = 0q. 
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Since ( 1.1 1 appears to be better behaved when /3 > 1, it is reasonable to expect 
that the previous theorem continue to hold in the entire range a, /3 € [1,2], cf. [U 
[T2J [131 HH] for some relevant results. 

In what follows we shall also need a local-in-time existence results for Sobolev 
regular solutions to the inviscid version of ( |1.1[ ). However, as mentioned in the 
introduction, the currently available results [3 [71 [2] apply only to the inviscid 
quasi-geostrophic equation {(3 = 1) 

6t + u-Ve = Q, u = (-7120,7110), 

with initial data 6*14=0 — 6*0 belonging to some Sobolev space H''. Throughout this 
paper we will simply make the standing assumption that the inviscid generalized 
quasi-geostrophic equation {6t + u- \79 = 0, u = curl A^'^ 6* for /3 G [1, 2]) possesses 
such a sufhciently regular solution. 



2.2. Fractional calculus. The fractional Laplace operator A" occurring in (1.1) 
is defined using Fourier transform, namely 

Our normalization of the Fourier transform reads 



JR2 Jr2 



In the upcoming analysis, we will need a different representation of A". Since we 
require a G (0,2), 7\ provides the identity 



|0|2+" ' " 7r2"r(a/2)' 

for all / in the Schwartz class (in particular all / G C^). Here, P.V.J ■ ■ ■ dz denotes 
the principal value integral, and we have introduced the notation 

We will make use of the following Leibniz-like formula "with remainder" : 
G"(/,g) :=A"(/.g)-/A"5-5A"/. 
By adding and subtracting, we see that 

Szifg)^gSJ + f6,g + SJS,g. (2.1) 



Multiplying (2.1) with \z\ ^ ^ and integrating over z provides the following repre- 
sentation of . 

Lemma 2.3 (Leibniz formula). // / and g are sufficiently smooth functions the 
following identity holds pointwise, 



G^{f,g)^K'^{fg)-fK"g-gK''f 



fC^P.V.^IP^S dz, /or a G (1,2), 
\2V/-V<?, for a ^2. 



In our analysis, we will make several applications of the above Leibniz rule. The 
following proposition provides an estimate for the term G". It is a variation of 
a result due to Constantin [3]. 

Proposition 2.4. Let 1 < p < oo. For each fixed a G (0, 2], there is a constant C 
depending on a such that 

||G"(/,.9)|Up <c(||v/|Uo.||v.9|U.)"/'(||/IUoo||5|u.)i-"/2, 

for all f, geC^. 
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Proof. The result follows directly from the Holder inequality when a = 2. We may 
thus assume that a G (0, 2). Let us write G"{f, g) as the sum of two parts: 



G"{f,g)=P.V. [ 



l\z\<r 

:=J^{f,g)+K^{f,g). 
(i) We commence by estimating the first term: 

p 

SzfSzg 



L 



SzfS^g 



\z\>r 



\2+a 



dz 



\Jrif,9)rL 



L" 



\2+a 



dz 



< 



\z\<r \Z\ 

\5zfSz 



dx 



\z\<r 



\2+a 



dz dx 



< 



< 



\qi\<r 



\qi\<r 



'||V/||z,oo||V5||lp dq,...dqp 



= 11 -I i{^^f%fdci....dci,dx. 

U^J\q,\<r J\qA<rtJ{ 1*1 

Since YTi=i — ^^^^ '^^'^ apply the generalized Holder inequality to obtain 
\\J?{f,g)\\%< [ f[\qi\-^-"\\S,JS,,g\\LP dq^...dqp 

J\qi\<r J\qp\<ri^i 

I i{\<li\-^-''¥qJ\\LAK9\\L^ dq,...dqp 

•'l9pl<'" i=l 

''\qp\<r- j=i 

\J\z\<r ) 
<,p(2-a)c||V/||i^||V5||L, 

where we have used that < |^|||V/||lp, for all / e W^^^. 

(ii) A calculation similar to the previous yields 

<c||/lliooNli. ( / 

\j\z\>r 

<r-''''c\mu\gr^,. 

Optimizing in r yields the result. 



< 



|9l|>r J\qp\>r 



-2-a| 



\lp dqi ...dqp 



-2-a 



dz 



□ 



In order to work efSciently in our Sobolev spaces we will need to provide some 
standard definitions, mostly to fix the notation. Let I denote a two-dimensional 
multi-index, i.e., / = {li,l2), Ij G No- Then we write 

0\i\f 



If ^ e N we let 



dx^^dy^^' 

vV = {vV I |Z| = 
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and 



9- 



\i\<i 



We will be working the Sobolev spaces 

iJ^- - //'^(M) = {/ e 5' I (1 + G (K)}, fc e M 

(where S' denotes the set of tempered distributions). If /c is a natural number, 
is the standard Sobolev space with inner product and norm given by 

k 

where we have introduced 

In our analysis, we will apply the following corollary of the previous proposition. 

Corollary 2.5. For each fixed a G (0, 2] and integer fc > 3, there is a constant C 
depending on a such that 

\\G'^{f,g)\\H''<C\\f\\H.+4g\\H.+., 

for all f, g e H''+\ 

Proof. Since C§° is dense in H^~^^, we may assume that f,g& Cg°. By definition, 



2 

L2- 



s=0 



Our strategy is to prove the desired estimate for each of terms in the sum sepa- 
rately. For this purpose, we let s = 0, . . . , fc be arbitrary and consider an arbitrary 
component of V^G°'{f,g). Using Lemma 2.3 and the standard Leibniz rule, we 
deduce 

-J^G"(/,.g) 



p.v. 



; s-i 



EE 



; s-i 

EE 



n— m— 
/ Qm+n 

\dx'^dy"' 
s - I 



s-l 
m 



dj 



s — m—n 



d 



(2.2) 



dz9 dz 



d 



\7n-\-n 



d 



s—m—n 



dx^dy"^ ' dx^^"-dy^~^ 



9 



n— m— 

For m, n such that 2 < m + n < /c, Proposition |2.4| can apphed to conclude 



G' 



a ( ^ 



d- 



■s—m—n 



\dx'^dy™- ' dx^~"dy'^~^ 



^9 



< C||/||^...||V™+\9||^/i||V~g||^-"/^ < C\\f\\H>^^^\\g\\H>^^^. 
Conversely, if m,n is such that < m + n < 1, Proposition |2.4| allow us deduce 



G 



( d 



d 



S—m—n 



L2 



<C||.g||H.+i||V' 



+Ji+1 j-||"/2||T7m+n r\il-a/2 



71 



/II 



< C'II.9|Ih''+i||/IIh* < C'II5||//'=+i||/IIh'--+i, 
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since fc > 3. Now, by taking the norm on both sides of (2.2) and applying the 
previous calculations, we gather 

QS 



dx^dy 



L2 



Hence, any given component of V^G°'{f,g) satisfies the desired bound. Thus, 



which completes our proof. 



□ 



2.3. Two auxiliary lemmas. We will make heavily use of the following two lem- 
mas throughout the paper. Their proofs are technical and tedious, but straightfor- 
ward. For this reason, proofs are deferred to the appendix. 



Lemma 2.6. Let k > 6 be an integer. Then, 

k 



for all f e H 



s=0 
k 



<C||/||H-HI/lPff., 



(2.3) 



Lemma 2.7. Let k > A be an integer. The following estimates hold 

|2 J- „ ^ rrfe 



V / V (V/ • curl A-'^.g) : V/ dx 
y] / (V.g • curl A-*^/) : V/ dx 



<C\\g\\H^+i\\f\\l„ feH',geH 



(2.4) 



fe-i-i 



(2.5) 



3. GODUNOV SPLITTING 



In this section we prove rigorously the expected linear rate of convergence for 
Godunov splitting. As part of the proof we also show that this splitting method 
is well-defined and that it produces regular approximations. The results are valid 
under a condition on the length of the time step At. 

We begin by precisely defining Godunov splitting for (1.1 1. For this purpose. 



write ( 1.1 1 as 



7* = 0^,7;, cie)^Ai0) + Bie), 

B{9) = -Ve, A{e) = -A'^O. 

Using the operators A and B, we define the solution operators $^ and $b as the 
solutions to 

dt^Ait, 60) - A($A(t, do)), $a(0, 60) = Bo, 

at$s(i,0o) -S($B(t,0o)), $s(O,0o) = 00. 

The Godunov splitting method is classically defined as follows: For At > given, 
construct a sequence {6*", 6'"+^/^}|j"^o'^*^ of approximate solutions to (1.1) by the 
following procedure: Let ~ 9q and determine inductively 

= $B(At,r), 61"+^ = $A(At,6i"+i/2), n = 0,...,[T/AtJ -1. 

To facilitate the convergence analysis, we will need a different definition of the 
Godunov method. Our definition can be seen as an extension of the splitting 
solution {6*", 0"^^/^}„ to all of [0, T]. The most used method of extension is to let 
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"time run twice as fast" in each of the sub-intervals [tn,t„+i/2] and [tn+i/2- ^n+i], 
where as usual = rAt for r G [0,c»), thus obtaining 



{t) = 



$B(2(<-t„),r), te [tn,ta+l/2), 



Although it appears to be a natural extension, it does not seem to be appropriate 
for our purpose. Instead, we will follow the approach taken in the paper jlOj . and 
introduce two time variables instead of one. 

Our Godunov splitting method is given by the following definition. 

Definition 3.1 (Godunov splitting). With A< > given, define the domain 

[T/AtJ-l 

^At = [J [^n,^n+l] X [tn,tn+l]- 

The time-continuus Godunov splitting solution § : f^At — > K is defined as the 
solution to 

Ut{t,t„) = B{d{t,tn)), t e (t„,t„+l], 

\^t=u = (3.1) 

Ur{t,T) = A{d{t,T)), {t,T) e [tn,tn+i] X (t„,t„+i], 

cf. Figure [T] 
Observe that 

^{tn,t„) = 0'\ n = o,...,Lr/Atj. 

Thus, i9(t, t) is an extension of {6'"}„ to aU of [0, T]. 
To measure the error, we will use the function 

where 9 is the (smooth) solution of 

It is not trivial trivial to obtain the existence of a splitting solution in the 
sense of Definition 13.11 Since the best available existence result we have for the 
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hyperbolic step is local-in-time, it is unclear if we can iterate the steps up to any 
given large time T. In fact, well-posedness of the method is one of our main results. 
The following theorem is our main result in this section. 

Theorem 3.2. Suppose 9q e with 6 < fc e N. and that a, jS € [1, 2]. Then, for 
At sufficiently small, we have the following: 

(1) The Godunov splitting approximation 'd(t,T) is well-defined. Moreover, 
d{t,T) belongs to C{[0,T]^; H''). 

(2) The error e(t) :— \\'i!}{t,t) — 9{t)\\ijk-2 satisfies 

||e(i)||^.-2 <tCAt\\eo\\jj,. 

Theorem |3.2| will be an outcome of the results proved in the subsections below 
(Subsection |3.3| will bring the pieces together). 



3.1. Evolution equations for the error. To prove Theorem |3.2[ we will analyze 
a set of evolution equations satisfied by the error e, see |10) . which we now derive. 
We shall need the following Taylor expansion satisfied by an operator E: 



1 



Eif + g)^E{f) + dEif)[g]+ / {1 ~ j)d'E{f + jg)[gY dj. 

Jo 

Using the definition of i9 and the above Taylor formula, we deduce 

et - dC{d)[e] ^^t+dr-0t- dA{d)[e] - dB{0)[e] 

= i9t + A{^) -{A + B){e) ~ dA{e)\e\ - dB(d)\e\ 

^dt- B{d) + {A{i)) - A{9) - dA{e)[e]) 

+ {B{d) - B{e) - dB{6)[e]) 
1 



(3.2) 



= F{t,t)+ / (l-7)d^C(e + 7e)[e]^ d7, 
Jo 

where we have introduced the "forcing" term 

F{t,T)^Mt,r)-B(d{t,T)). 

By direct calculation, 

Fr - dA{d)[F] = Vtr ~ Bid)r " dA{'d)['&t " B{d)] 

= A{d)t - dBid)[^r] - dAi^)['dt] + dAi^)[Bid)] 

^ dA{d)['&t] - dB{'d)[A{i))] ~ dA{i))[^t] + dA{'d)[B{'d)] 

^[A,Bm, 

where we have defined the commutator 

[A,B]{f) = dAif)[Bif)] - dBif)[Aif)]. 



(3.3) 



For the fluid equation (1.1), we have that 

Aif) - -A"/, 
dA{f)[g]^~A^g, 
d^A{f)[g,h] = 0, 



and 



S(/) = -V/-curlA-'^/, 
dB{f )[g] = -V/ • curl A"^.g - Vg ■ curl A"'^/, 
d^B{f)[g, h] = -Vh ■ curl A~^.g - Vg ■ curlA^^h. 
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The Leibniz formula (Lemma |2.3[ ) yields 

[AS](/)=A" (V/-curlA-^/) 

- V/ • curl A-'' A"/ - (VA"/) • curl A^'^/ 
= G"(V/,curlA-^/). 



Hence, for the fluid equation (1.1), equations J3.2|-(3.3) read 



F = i?f +Vi9-curlA-''i?, 
Fr+A°'F = G"(Vi?,curlA-'^i?), 
et+A^e^F-Ve- curl A"'' 6* - V9 ■ curl A^'^e - Ve • curl A~^e. 



(3.4) 
(3.5) 



The equations (3.4 1 and (3.5) constitute our main tool for proving Theorem 3.2 



3.2. Estimates on the error. Lemma 3.5 below provides an estimate on the error 



e that will be a key ingredient in the proof of Theorem 3.2 The estimate depends 



on two auxiliary results (Lemmas 3.3 and 3.4 1 that we first prove. 

The first auxiliary results gives perhaps the most fundamental property of our 
splitting solution. It states that if the splitting solution is in H'', fc > 6, then it is 
actually in H''+^. 

For notational convenience, let fi^j denote the set of all times prior to {t,T): 
^*At = {(s^^) e ^^At I < s < i, 0<<j<t}^ f^At n ([0,t] X [0,r]). 
Lemma 3.3. Assume that a, /3 G [1,2], Q < k E N, and let 'd be the Godunov 



approximation to (1.1) in the sense of Definition 3.1 and (3.1) 



If for some {s,a) G 57a, 



At' 



(3.6) 



then 



Proof. By direct calculation we see that 



drlMt,T)\\l, = V / y'^r ■■ V'?? dx 

k 



- 1 



+ a/2- 



Hence, for (i,T) e [t„,t„+i] x (i„,t„+i], n = 0, . . . , IT/At\ 

mt,T)\\l, < mt,t,,)\\l,. 
We now calculate a bound on the H'^ norm of i?(<, i„). By definition, 

= - V / (Vi? • curl A-'^ J?) : V^d dx. 



(3.7) 



Applying Lemma |2.6| to this identity, we see that 

dt\mt.t^)fH. < cmt,t^)\\H>^-4d{t,t^)\\i,. 

Clearly, this allows us to conclude that 

dtmt,t„)\\H^ <cj\\{){t,t„)\\H^, 



(3.8) 
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and integration in time gives 

Consequently, 

which completes the proof. 

The following lemma is our second auxiliary result. 
Lemma 3.4. Let k, t, and r he as in the previous lemma. Then, 



\F\\%k-2 + a/2 



□ 



Proof. Applying V" to ( |3.4| , multiplying componentwise with \/^F, summing over 
s = 0, . . . , fc — 2, and integrating over the domain gives 

fc-2 

(3.9) 



= V / : V'G"(V79,curlA-''z9) dx 

< C||i^||H'=-2||G"(Vi?,curlA-'3i9)||^.-2. 



An application of Corollary |2.5| to (3.9) gives (here we do not make use of the 
regularizing effect of A~^) 

drl\\F\\l,-2 < C\\F\\H^-4S/nH^-4cmlA-"nH>'-^ < C||F|1^.-. ||^||^.. 



In view of Lemma |3.3[ this means that 
Since F{t,tn) = integration in time gives 
which concludes the proof. 



□ 



The next lemma will be the key ingredient in the proof of Theorem 3.2 
Lemma 3.5. Let k, t, and t be as in Lemma \3.3\ Then, 

||e(s)||fffe-2 < sAiCie='^*||6lo||^fc, < s < t. 



Proof. Applying to (3.5), multiplying componentwise with Ve, summing over 
2 



s — 0, . . . ,k ^ 2, and integrating gives 
1 
2 



fc-2 



s=0 



= J2 V^F : V"e + V" (Ve-curlA^'^e) : V' 



V (Ve* • curl A-'^e + Ve • curl A-'^6l) : Ve 

fc-2 



dx 



(3.10) 



< e 



Hk-2\\F\\H>^-2 +Y] f I V" (Ve-curlA-'^e) : V"e 

fc-2 

- J2 (^^ • curl A-'^e + Ve • curl A^^^ 6*) : V"e | dx. 



s=0 
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An application of Lemma 2.6 provides the estimate 

k-2 



V / V (Ve • curlA-'^e) : dx < C\\e\\]j,^2 



s=0 



(3.11) 



where we have also used that ||i?||//fc-2 < 7. 

By virtue of Lemma |2.7[ we also have the estimate 

fc-2 



V / V" {ve ■ curl A-^e + Ve • curl A-^O) -.V'edx 



(3.12) 



< 



c {\\e\\H>^-^\\e\\l,^. + \\e\\H.-A\4]i^-A 



Combining (3.12), (3.11), and applying the result together with Lemma 3.4 



to 



(3.10) enable us to conclude that 

An application of the Gronwall inequality (recalling that e(0) — 0) yields the result. 

□ 



3.3. Proof of Theorem |3.2[ We will make use of the following bootstrap lemma 
(cf. [211 Proposition 1.21]): 

Lemma 3.6. For each (t,T) G ^At, suppose that we have two statements, a "hy- 
pothesis" H{t,T) and a "conclusion" C{t,T). Suppose that we can verify the fol- 
lowing assertions: 

(1) If H(t,T) is true for some {t,T) € ^At, then C{t,T) is also true. 

(2) If C{t,T) is true for some {t,T) G ^At, then H{t' ,t') is also true for all 
{t' ,t') in a neighborhood of(t,T). 

(3) //(ti,Ti), (^2,7-2), ••• is a sequence in Q,At converging to {t,T) £ flAt, and 
C{tn,T„) is true for all n, then C(t,T) is true. 

(4) H{t,T) is true for at least one (t, t) G f^At- 
Then C{t, r) is true for all (t, r) G ^At- 

Let H{t, t) denote the statement 

ll^(s,0-)||//fc-2 < 7, (s,cr)erJ^t, 

and C{t, t) the statement 

\\^is,cr)\\H>,-2 (s,cr)er2^;, 

where 7 is some value which will specified below. Let us for a moment assume that 
assertions (1)~(4) of Lemma 3.6 are true for H and C. In this case Lemma 3.6 tells 
us that C{t,T) holds for all {t,T) g ^At- Hence, the H''~^ norm of the sphtting 
solution never blows up and thus the existence part of Theorem 3^ follows readily 
(a local-in-time solution can be extended to all of (0,T)). 



Let us now verify assertions (l)-(4) of Lemma [376| for our choice of H(t, r) and 
C{t,T). First, we observe that assertions (2) and (3) clearly hold. For assertion (4) 
to be true, we assume that 7 satisfies 

\\mo)\\m-2 = \\9o\\m-2 < 7- (3.13) 

Then it remains to verify assertion (1). For this purpose, let us assume that H{t, t) 
is true for some (t, r) £ i^At- Lemma |3.5| can then be applied to obtain the bound 

\\e{s)\\H^-2 < s{At)Ce'^'\\eo\\l,, s e [0,t]. (3.14) 
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We need to compare i?(s, cr) with the corresponding value on the diagonal, viz. 
i9(s, s). Two cases need to be considered: 
(i) If cr > s, we observe that 



s=0 • 



s=0 • 



(using that a < 2 to make sure the expression is finite) which implies that 

ms,a)\\H.-2 < Ms,s)\\h^-2. 
(ii) If (7 < s we first see that 

fc-2 „ 



s=0 
fc-2 



V / \/'A°'d : V'l? dx 



< ll'?(S:0-)llH'=-2||'9(s,(T)||^fc, 



(using a < 2) which implies that 

Using this inequality and applying Lemma |3.3| we find 



\\^{s,a)\\H,-2 <\\dis,s)\\H>=-2 + / \d4^{s,a)\\\da 

J cr 

< ||i?(s,s)||^fc-2 + At sup \\d{s,a')\\Hi^ 

o-'e [ct.s] 

< ||z?(s,s)|U.-.+Ate^^^||0o|U.. 
Thus in both cases we conclude that 

||t9(s,a)||Hk-2 < ||i?(s,s)||^fe-2 + Ate'^''"||6io||//fc, |s - cr| < At. 



By applying the previous inequality, adding and subtracting 9, and involving (3.14 1, 
we estimate 

||^(s,a)||^.-2 < \\eis)\\H.-2 + \\9{s)\\h.-2 + Ate^-''\\0o\\H^ 

< CAte^^^(s||0o|lH. + \\0o\\h^) + C\\0o\\h^ 

< AtCi(7) + C2. 
Now, we fix 7 and At according to 

C2 



(3.15) 



7 = 4C2 



At < 



Ci(7)' 



and note that this is not conflict with (3.13). 



Then, (3.151 gives 



||z?(s,a)|U.-2<2C2 = ^, (s,a)ef]^;. 



which verifies (1) in Lemma 3.6 



At this point we have verified assertions (l)-(4) of Lemma 3.6 for our choice of 
H{t,T) and C{t,T). Consequently, Lemma [3. 6| tells us that C{t,T) is true for all 
times (t,r) G f^At- In other words. 



which concludes the existence part of Theorem |3.2| 



(3.16) 



Equipped with (3.161, we can apply Lemma 3.5 to obtain the error estimate 
||e(t)||H.-2 <tAtCe^*||0ol!?^., [0,T], 



which is the second part of Theorem 3.2 



□ 
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Figure 2. A schematic view of Strang splitting. 



4. Strang splitting 
In this section we prove the expected second-order convergence rate of Strang 



sphtting appUed to (1.1). We also prove that the method produces regular solutions 



up to any given finite final time T > 0. The results are valid under a condition on 
the size of each time step At. 

We will continue to use the notation and definitions introduced in the previous 
section, unless explicitly stating otherwise. The Strang method we consider in 
this section reads as follows: For At > given, construct a sequence {6'"}|j^o'^*^ 
of approximate solutions to (1.1) by the following procedure: Let 6*" = 9q and 
determine sequentially 



nn+l 



As with the Godunov splitting, the convergence analysis will require a time- 
continuous interpolation of the Strang approximations {0"}n- In contrast to the 
Godunov case, we will now introduce three time variables instead of two. This 
approach is different from the one taken in jlOl , and indeed appears more natural. 
In [lOj the authors stick to two time variables and interpret Strang splitting in 
terms of two "Ai/2" Godunov splittings, and in alternating order. 

Definition 4.1 (Strang splitting). With At > given, define the domain 



LT/AtJ- 
ri=0 



til t 



ri+1 



+ 1J 
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The time-continuous Strang sphtting solution d : f^At 
t9(0,0,0) = ^?o, 

'dr{t,T,'-f)^A{d{t,T,'-f)), (i,T)e 
'l?^(t,T,C^) =B(l?(t,T,L^)), {t,T,Oj)£ 



is defined as 



2 ' 2 



2 ' 2 



X (in ; ^ri+l] 7 



(4.1) 



cf. Figure [2| 
In each box 



2 ' 2 



n+lj 



2 ' 2 



, it is of particular interest to 



consider the function z9 restricted to the diagonal, i.e., the function ^ |) for 
t € [tn,tn+i]- Observe that at each point on the diagonal, is a Strang splitting 
solution with a specific time step. Specifically, is a Strang splitting 

solution with time step t — tn- Consequently, 



= 0,...,LT/AiJ, 



and hence that i? (|, i, |) can be seen as an extension of {6'"}„ to all of [0, T]. 

To measure the error between the splitting approximation and the exact solution, 
we will use the function 



where 9 is the (smooth) solution of ( |1.1[ ). 
We state two main results in this section. 



0(i), 



Lemma 4.2 (Well-defined). Let -d — i9(t,T,uj) be the Strang splitting solution of 
(1.1) in the sense of Definition 4-1 Suppose 9q € with 6 < fc € N, a, /3 € [1,2], 
and that At > Q is sufficiently small. Then 

(1) the Strang splitting method is well-defined with d G C([0, T]^; i?'^); 

(2) the error e{t) — i9(|,t, |) — 0{t) satisfies 

\\e{t)\\Hk-2 < CtAt. 
for some constant C — C (||^ollff'«)- 
Theorem 4.3 (Convergence). Under the conditions of the previous lemma, 

||e(i)||^.-3„ <C(Ai)2, 
for some constant C = C (||0oll^r'=)• 



Lemma|4?2]and Theorem |4.3| will be consequences of the results stated and proved 
in the ensuing subsections. 



4.1. The error evolution equations. To prove Lemma [4~2| we will use the same 
approach as for the Godunov method. However, since we are now using three time 
variables instead of two, we must derive a new set of evolution equations governing 
the error. 
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By a direct calculation, the error e satisfies the time-evolution 
et - dC{e)[e] = '^+dr + '^~et~ dA{e)[e\ - dB{e)[e] 

= y + + \B{d) -{A + B){e) - dAi9)[e] - di?(0)[e] 

+ (Aid) - A{0) - dA{9)[e]) + {B{^) - B{e) - dB{9)[e]) 

^F{t)+ f (1 - 7)^2^(0 + 7e)[e]2 d7, 
Jo 

where F{t) = f ) and 

F{t, T,Lu)^^ (Mt, r, uj) - B{-d{t, r, uj))) + - A{d{t, r, uj)). 

Since dr — A{d{t, r, uj)) — when lo — n = 0, . . . , [T/AtJ — 1, we can apply 
the arguments of (|3.3|) to obtain 



F,-dA{d)[F] = ]^[A,B]{d), ef^At, 

where 

[A,B]{f ) = - dBU)[AU)]. 

We also derive the following equation for the evolution of in w: 

F^-dB{d)[F] 

= \^t. - \B(i))^ - \dBmdt - B{d)\ 
+ I?™ - - dB{§)\dr - A{^)\ 

= \B{^)t - \dBm^u] - \dBmM + UBmBm 

+ dB{i})[dr] - dA{d)[B{^'))] - dB{d)[dr] + dB{i})[A{d)] 

= ^ {dB{d)[§t] - dB{d)[B{d)] - dB{d)[dt] + dB{§)[B{§)]) 

+ dB{d)[A{d)] - dA{d)[B{d)] 
= dB(§)[A{d)] - dA{d)[B{d)]. 

Thus, recalling the definition of [ • , • ], 

F^^dB{^)[F] - [B.Am. 



For (1.1 1 the above evolution equations take the form: 

= ^ + Vi? -curlA-^??) + A"i9, (i,r,a;) e rjAt, 

+ A"^' = iG"(Vi9,curlA^^7?), a; = ^ , n = 0, . . . , [T/AtJ , (4.2) 
= -Vz9 • curl A^^F - VF • curl A"'' 

- G"(V79,curlA~'^i?), (i, r,a;) e l^At, (4.3) 
et + A"e = F - Ve • curl A^'^g - • curl A~'^e 

- Ve • curlA^'^e, (t, r, w) e f^At. (4.4) 



To prove Lemma 4.2 we will adapt Lemmas 3.3 3.4 and 3.5 to the Strang 
splitting approximation. 
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Let f^^j''^ denote the set of times prior to (t, r, w): 

■\t,T.UJ 



n 



At 



^Atn{[0,t] X [0,t] X [0,uj]). 



(4.5) 



We begin by observing that Lemma |3.3| can be easily adapted to yield the fol- 
lowing result. 

Lemma 4.4. Let 6 < fc G N, and assume the existence of [t, r, uj) £ f^At such that 

(4.7) 



Then 



||79(s,a,C)||H'^-2 <7, (s,a,C)ef^L ■ 



Next we prove a Strang version of Lemma 3.4 



t.T,UJ 

At ■ 



Lemma 4.5. Under the conditions of the previous lemma, 

\\F{s,a,C)\\H^-2 < AtCie^=^(^+«), {s,a,0 € ^ 
Proof. First, we observe that in the plane given by cj = 

F{s, i^t - , n^O,..., LT/AtJ - 1. 

Note that while ^ on the line (s,i„, ^), there holdsJ_(i?t - B{-d)) = on 



{s,tn, -^). Hence, we can repeat the arguments of Lemma 



3.4 



to conclude that 



||F(s, a, ^)||^.-2 = II- l^M^, a, p - i?(^(s,a, ^)) ) ||^.- 



(4.8) 



We fix n such that {t,T,uj) e [^, ^^^] x [i„,i„+i] x ^^ii]. Now, to estimate 
||i^||//fc-2 at an arbitrary point (s,cr, C) G ^At "' ^® ^^^^ ^'"^^ integrate in the uj 
direction to the plane given by w = ^ and then apply the estimate (4.8 1. To 
perform this integration, we apply V* to (4.3 1, multiply the result with VF, sum 
over s = 0, . . . , fc — 2, and integrate by parts, to obtain 

fc-2 



3c.;^||i^||^.-2 ^-Y [ V^F : V^G"(V79,curlA-^i?) dx 

fe-2 

y / V" (Vi? • curlA-'^i^) : V"F dx 
y / V" (Vi^ • curl A-'' t^) : da; 



(4.9) 



s=0 ■ 

:=/i+/2+/3, 

and the goal is to bound the terms /i, I2, and /a. 

By first estimating as in (3.9), then using Corollary 2.5 and finally applying 
Lemma 4.4 we obtain 

|/i| < ||F||^.-2||G"(Vi9,curlA-''i9)||H.-2 

< C||F||afc-2||Vi?||jjfc-i||VcurlA-'^i?||Hfc-i 

<Gi||^^||^.-2e^^^(^+?)||0o|lH.. 

To bound the two other terms I2 and /a, we first apply Lemma |2.7| to find 

I/2 + /3I <q|F||^._2||i9||H.-i. 
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An application of Lemma 4.4 then yields 
Inserting these bounds back into ( |4.9[ ), 
Applying the Gronwall lemma we obtain 



\Fis,aX)\\H''-^ < e 



Hk-2 



Finally, we apply (4.8 1 to conclude 

\\Fis,aX)\\m-^<^tC\e^'^'<^'+'^\ 
which completes the proof. 

Lemma 4.6. Under the conditions of Lemma \4-.4\ 

\\e{s)\\Hi'-2 < sAtC{j), se[0,i], 
where C{-y) is a constant depending only on j, T, and ||^o||fffc- 



□ 



Proof. Since (4.4) is identical to (3.5), we can repeat the proof of Lemma 3.5 step 
by step (using Lemmas 4.4 and 4.5 instead of Lemmas 3.3 and 3.4 1 to conclude. □ 



4.2. Proof of Lemma 



4.2 



Recall that ^l^^j'" denotes the set of all times prior 
to (i, r, oj), cf. (4.5 1. To prove the existence part of Lemma [4. 2| we will utilize the 
same strategy as with did for the Godunov method in the previous section. That 
is, to apply the bootstrap lemma (Lemma 3.6) with H{t,T,uj) as the statement 

ms, a, OWh"-^ < 7, (s, a, Q e f!^;'", (4.10) 
and C(i, t, uj) the statement 



Ms,aX)\\m-2 < 



7 



where 7 is some value to be specified below. Note that there is no problem with 
extending Lemma 3.6 so that it allows three time variables instead of two (i.e., 
hypothesis and statement depending on three variables). 

To conclude the existence part of Lemma 4.2 we must verify assertions (l)-(4) in 
Lemma 3.6 Assertions (2) and (3) clearly hold. Regarding assertion (4), we make 
the assumption 

11^9(0,0,0)11^.-2 = |10oI1h^-2 <7- (4.11) 
It remains to verify assertion (1). For this purpose, let us assume that H(t,T,Lu) 
is true for some (i,T, w) e ilAt- In order to show assertion (1) we need to prove 



that C{t,T,uj) is true. We first note that H{t,T,uj) being true renders Lemma 4.6 
applicable, allowing us to conclude 



le(s)|l^fc-2 < sAiCi(7), s e [0,t] 



(4.12) 



Now, to estimate ||i9(s, ct, C)|lHfc-2 at any time (s,(t, C) € f^^j''^ the idea is to 
consider the time evolution of 1119(5, a, C)||^fc-2 in the oj and r directions. Let n be 



such that (t, T, w) G 



First we integrate from (s, a, C) 



to (s,cr, -f-). Then, we integrate the result in the plane w 



from (s, cr, ^) to 



(s,2s, ^). Finally, we integrate in the uj direction from (s,2s, ^) to the diagonal 
(s, 2s,s). On the diagonal, we will then utilize (4.12) to conclude the necessary 
bound. 
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By definition, 



1, 



s=Q 
k-2 



s=0 



^ / V {Vi9 ■ curl A-^i9) : Vi? dx 



(4.13) 



where the estimate follows from Lem ma |2.6[ and where we have used (4.10). The 
fundamental theorem of calculus and (4.131 provides us with the estimate 



(s,cr,C)||H'=-2 < lll?(s,cr,C2)ll_ff'=-2 + CAt sup ||l9(s,cr,C')llH'»-2 



(4.14) 



<ms,a,C2)\\H>'^-+CAt\\9o\\H^-2, 



for any 'f < C2 < 
Similarly, 



9. 



1 



d [ s,a, 



k-2 



s=Q 
k-2 



J2 I "^'^T ■■ V"?9 dx 

■0 

J2 I V^A^i? : Vi? dx 



< 



s=0 



where we have used that a < 2. Thus 



d-. 



1 



tr. 



JJk-2 



< 



Hk-2 



d \ s,a, 



d [ S,(T, 



From the fundamental theorem of calculus, using the previous inequality and Lemma 
14.41 we conclude 



tr. 



< 



Hk-2 



s,2s. 



< 



Hk~2 

At sup 

ct'G [min{ 2s, cr} ,max{2s,cr}] 
tn 



(4.15) 



d s,2s, 



Thus, by combining (|4.14|) and (|4.15|), 

i9 



(s,c^,C)IIh'=-2 < 



t. 



< 



Hk^2 



Uk 



d s,2s, 



-A^CIl^olU- 

CAt\\0o\\H^-. 
f AtC||^?o||ff'= 



Hk-2 

< \\^is,2s,s)\\j,k-2+AtC\MH>^ 
<Ci(7)Ai + C2, 

where the last inequality is ( |4.12[ ). Now, we fix 7 and At according to 

C2 



(4.16) 



7 = 4C2, At< 
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and note that this is not conflict with (4.11). Then, (4.16) gives 



7 



3.6 



We 



Since {s,a,(^) € ^At' arbitrarily chosen, this verifies (1) in Lemma 

have now verified assertions (l)-(4) of Lemma 3.6 for H{t,T,uj) and C{t, r, lu) 
Consequently, Lemma 3.6 tells us that C{t, r, lu) is true for all times {t, r, w) £ ilAt 
In other words, 



\\d{t,T,io)\\H>^-2 <j, {t,T,uj)enAt- 



(4.17) 



Thus, we can conclude the existence part of Lemma 4.2 Since (4.17) holds, we can 
apply Lemma 3.5 to conclude the error estimate 

\\e{t)\\Hi^-2 <tAtC, te[0,T], 

which is the second part of Lemma |4.2| 

□ 



4.3. Temporal regularity and auxiliary estimates. So far we have focused 
on the spatial regularity of the splitting solution. To prove Theorem |4.3[ we will 
need regularity in time of z9. Since we have introduced three time variables, it is 
convenient to define a time gradient. We will use the notation 



ftt 


ftr 


ftUJ \ 


frt 


frr 


fruj 1 


fujt 


fujT 


fuJUJ / 



with the obvious extension to higher order (i.e., / for fcth derivative). 

Lemma 4.7. Let a d [lj2] and be the Strang splitting solution in the sense of 
Definition 4-1 Then, for any k and I such that 9o & 



Proof. We argue by induction on /. For / = 0, the result is merely an application 
of Lemma [4. 4| Now, let us assume that the bound holds for ^ = 0, . . . , g. To close 
the induction argument, it remains to show the bound for I — q + 1. 
Let (t',T',uj') € ^At be arbitrary and fix n such that 



{t',r',oj')e 



^ [^n: ^n+l] ^ 



2 2 

An arbitrary component of V^^^i) can be written in the form 

Qq+l 



0<i,j,£<9 + l, i + j +£ = q + 1, i,j,i&N. 



To estimate the arbitrary component at the point {t',T',uj'), we first apply the 
fundamental theorem of calculus to write 



1, 



2 



k 

E 

s=0 • 



r I V^a„e?+i(i',T',s) : V^e?+^,(i',r',s) dxds. 
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By applying the definition of the definition t?^ (cf. (4.1)), and tlie Leibniz 

rule, we deduce 



/ ^'d^e'+]{t',T\s) : ^'e'+]{t',T',s) dxds 

JR2 



s=0 ~ 
k i j I 

s— r— n— rn— 



59+1 



r J \ n I \m 



: V''e«+^^(t',r',s) dxds. 
Inserting this expression into the previous equality yields 



1 

_ 1 
^ 2 



y -.^ ■. 2 

-EEEE 

s^O r— n— m— 



(4.18) 



/ \n J \iTi 



(curiA-^e;:;„+;r(t',r',.) . vet;-";"_™(t',r',5)) 

: V'0fr^(*'-^'.^*) dxds. 

Let us consider three separate cases of m + n + r in the quadruple sum above, 
(i) If TO + n + r = (7 + 1, the corresponding term in the above reads 

E r / (curlA-'3e,?+j(t',r',5) • Vi9(t',T',5)) : V^e?+;(i', r', s) dxds 



<C / ||^(i',r',s)||^...||e^+i,(t',r',s)|l^. ds 



(4.19) 



Here, the first inequality is an application of Lemma |2.7| and the last inequality is 
an application of Lemma |4.4[ 



(ii) If TO + n + r — 0, we can also apply Lemma 2.7 to conclude 

E r /, (curlA-^i9(t',r',s) • Ve,?+}(t', r', s)) : V'Q^+^it' ,t' , s) dxds 



<C / ||^(i',r',s)|U.+i|ie?+i,(t',r',s)f^. ds 



<c\\eo\\H^.^ / l|e^+;,(t',T',s)||^ 



ds. 



(iii) For the remaining cases (1 < m + n + r < q), we first apply the Leibniz rule to 
the spatial derivatives, then we overestimate the terms using the Holder inequality. 
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In this way, we obtain 

/ v^(curlA-%;^'„+;;+'■(^^r^.).ve?+l-:;— 

: V'e«+](i',r',s) dxds 

lie?+;,(t',T',s)||^.||e™+:;+'^(t',r',s)||^.., (4.20) 

<cj^^ \\e1+^,{t',r',s)\\H4eo\\l.+.i,+r> ds, (4.21) 

where we have used the induction hypothesis and a > 1 to derive the last inequahty. 
By applying (|4.19|)-(|4.21|) to (|4.18|, we gather 



i||e?+^,(t',r',..')ll 



2 



< 



+ C||0o||^.+(i+„. i m^!i{t\r',s)\\Hk ds 



+ C\\eo\\H^^^ / m^}{t',r',s)\\l, ds 



< 



09+1 ( t' r' — 



Hk 



+ CAt + C I ||e^'+}(t',T',s)|j^. ds 



where the constant C depends on ||6'o||//fc+(i+<i)Q ■ By applying Gronwall's inequality 
to the previous inequality, we conclude 



||e^+i,(t',r',c.')|ll,. < 



em t 



i I 4-1 ^" 



CM 



Hk 



(4.22) 



We have now derived a bound on 0^^^(t', r', w') in terms of Qj'^\{t' , t' , ^) . 

Next, we derive a bound on 9^^^(t',T', ^) in terms of 0^^^(t',t„, ^). For this 
purpose, we once more apply the fundamental theorem of calculus to obtain 



09+1 it' t' — 

_ 1 
~ 2 



^ £ v^9.e?+i, (t', s, I) : v^e,?+i, (i', s, ^) d^ds 



A; 

-E 



t,. JR2 2 \ 



2 



dxds 



ds, 
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where we have also used that {fr — A{'d) for w 

2 



09+1 ( ^' ^' h}L 



' ' 2 



< 



if-' 



It follows that 

2 



09+1 ^' / ^ 



(4.23) 



Finally, we perform our last application of the fundamental theorem to obtain 



1 



p.q+1 ( ^ 4- ^ 
V 2 2 



H'' 



By applying the same calculations to the previous equations as those used to 



derive (|4.22p, we arrive at 
t 



< 



\ 2 



1 / 

7 ''n 



(4.24) 



Combining (4.22|, (4.23), and (4.24) gives 



< 



e 



q+l I tn 



2 / _f^fc 



which immediately leads to 



N 



<||e?+i,(0,0,0)||^.e^*'+Ct'e"'. 



(4.25) 



Finally, let us estimate ||9^^^(0, 0, 0)||/jfc. By definition, we have that 



a«4 



At this point, we can apply each of the time-derivatives to d{0, 0, 0) and use the 
definition (4.1) to translate time-derivatives into spatial derivatives. Since a > 1, 
it is clear that the case j = q + l contains the highest number of spatial derivatives. 
In this case, 

d^+^d = -A"99?9 = (-1)«+ia(«+i)"^i. 
Thus, for any valid combination of and £, we conclude 



0?ti(O,O,O)||^. <C||0o||ff.+=(,+i). 



From this, and (4.25), we conclude that the lemma holds also for I = q + 1. 



□ 



The following lemma is almost a corollary of the previous lemma. 

Lemma 4.8. Let be the Strang splitting solution in the sense of Definition \4-T\ 
and (4.1). Then, for any k such that \\9Q\\fjk+3a < C, we have 

9o||/i-fc+3a) , (t, T,iLj) e JIa*- 



\VlF\\Hk <C\ 



Proof. Let i and j denote any one of t, r, or uj. An arbitrary component of Vj-F 
can then be written Fij :— didjF . By definition, we have that 

'1 



Fij = didj 



{-dt + Vi9 • curl K-'^-d) +dr+ A«^? 



= \dt'd,j + dr^ij + i cm\K-f^Vdij -Vd +]: curl A^'^Vi? • Vi?, 



^ curl A-^Vz9j ■\/^^ + ]^ curl A"'^V?9i • + A^'d. 



I] ■ 
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By applying the Holder inequality together with the previous lemma, we estimate 



+ i||V79, • curl A-'3i9,||^. + i||V79, • curl A-'^i^J^, 
< C+ ||V'''+i|?||L=e||V?79||H'= + ||V?i9||h'=||V'=z?|1loo 

+ 2||Vtt?||^.||A-'3Vti9||H^+3 
<C{l + \\d\\ 

where we h ave grossly overestimated most of the terms. We have also applied 
Lemma [4.7l with Z = 1 for the H^+'^ norm and with / = 2 for the norm. The 
constant C depends on ||6'o||Hfc+3a. □ 



Theorem 4.3 is a consequence of the following (remarkable) fact: 
Lemma 4.9. There holds 




Proof. We will prove Lemma |4.9| by direct calculation. Let us begin by estimating 
Fu (^,i„, ^). Since F{t,tnjf)^ 0' (EH* ^^l^s us that 



F I — f — 



-G"(Vi?,curlA-''79). 



Similarly, we see that (4.2) yields 



1 



Fr = -G"(Vi9,curlA-'3z9). 



(4.26) 



(4.27) 



It remains to estimate Ft (^,<„, ^). However, as F{t,tn, ^) = 0, we must have 
Ft{^,tn,*-f) =0. This, together with ( |4.26[ ) and ( |4.27[ ), concludes the proof. □ 

Using the previous lemma, we can now prove that the error produced along the 
diagonal (t/2,t, i/2) is second order in At. 



and (4.1). Then, for any k such that ||0o||_ff''+3a < C, 

< CiAtf. 



Lemma 4.10. Let be the Strang splitting solution in the sense of Definition \4.i\ 

(4.28) 



, t t 

F -,t, - 
' 2' ' 2 



Proof. Since i^(^,t„,^)=0, a Taylor expansion provides the identity 



, t t 

^'2'^'2 




ds (4.29) 
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where the last equality is an application of Lemma 4.9 By taking the H norm on 



both sides of (4.29) and applying the previous lemma, we gather 



, t t 



which concludes the proof. 



□ 



4.4. Proof of Theorem 4.3 We have now gathered all the ingredients needed to 
prove the sought after second-order error estimate. 



Performing the same calculations as in (3.10), (3.11), (3.12), and then applying 
Lemma [4.10[ yields 

i9t||e(t)||H^-3. < C{T) (At2 + \\e\\H.-..) , t e [0,r]. 

Since e(0) = 0, an application of the Gronwall inequality to the previous inequality 
gives 

||e(t)|Ufc-3c < iAt^C, ie[0,T], 
which concludes the proof of Theorem |4.3[ □ 



Appendix A. Proof of Lemmas 12.61 and 12.71 



The purpose of this appendix is to provide proofs of Lemmas 2.6 and 2.7 Both 



lemmas have been crucial to our convergence analysis. In particular. Lemmas 3.3 
[331 1431 and lO all rely on their validity. 



Lemma A.l. Let k > 6 be an integer. Then 

k 



V / V"(V/ • curl A-''/) : V"/ dx 



<c'|1/I1h-^I1./I1h., 



(A.1) 



for all f e H'' 



Proof. Let us confine to the three-dimensional case {N = 3) as the other cases are 
almost identical. 

By applying the Leibniz rule (with multiindcx notation a = (ai, 0:2, as)), we 
obtain the following expression 



V / V"(V/-curlA-'3/)V"/ 



\a\=s 



ai a.2 ctz 



EEEE 

\a\=s ■ji=0 12=0 43=0 

1^11-1-42+13 J 



Oi2 \ I 013 
h ) V «3 



■ curl A ^ 



d'f 



dx°'^dy°'^dz°''^ 



(A.2) 



dx. 



Let us now consider four separate cases of ii + 12 + *3 in the above quadruple sum. 

(i) If ii -I- 12 -I- 13 — s (i.e (01,02,03) = (ii,i2,i3) ), the above term can be 
rewritten as follows 

' i f^. • curl A-'^ f) f dx 



1. 



d'f 



dx°'idy'^^dz°'3 



■ curlA"*^/ dx = 0. 
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(ii) E 2 < ii + i2 + is < k - 3 

flil+i2+i3 f 



< 



Qii+i2+i3 j 



■ curl A 



d 



S — il— ■i2— ''■3 



/ 



dx 



curl A 



_I3 ( 5'*-H-»2-*3j 



dx^^dy'-^dz"-^ 

<C\\f\\H4f\\H>'-4f\\H>'- 

(lit) If fc < ii + 12 + is < fc - 1 



L2 



dx^^dy^^dz" 



L2 



< 



dx^^dy^^dz^s 

Qil+i2+i3 f 
1 " 

dx^^dy'^dz^^ 



■ curl A 



S — l]_—l2—'l3 



f 



g^ai-Hgya2-i2Qza3-i3 J J dx°'^dy°'^dz°' 



dx 



L2 



curl A 



Qs-ii-i2-i3 J 
52,ai-ii5yQ2-i29^a3-J3 



L2 



<^^ll/llif'^ll/lk^ll/llH'= <C||/||h-HI/IIh.. 

(iv) If < ii + ^2 + is < 1 

Qs-ii-i2-i3 J 



. V /. ^ ■ /■ -curlA-^ 
RN V dx^^dy^^dz^-^ 



9a;"i-'i9y"2-J252"3-*3 / / dx'^^dy'^^dz°'=> 



dx 



< 



gil+i2+i3 J 



dx^^dy^^dz^^ 



curlA-'^ 



Qs-ii~i2-i3 j 



^2,Ql-il^yQ2-i2 5z"3-«3 



L2 



5^/ 



dx^^dy^^dz" 



L2 



Hence, by applying (i)-(iv) in ( |A.2[ ), we see that 

V / V"(V/-curlA-'5/)V"/da:<C||/||ff.-2||/||^., 

which concludes the proof. 

Lemma A. 2. Let k > 4 be an integer. The following estimates hold 

k I. 

<C\\g\\HA\f\\]i>^, f,9^H\ 



V / (V/ • curl A-'' .g) : V/ dx 
y] / V'* (Vg • curl A-''/) : W dx 



(A.3) 
□ 

(A.4) 



<^^II5|Ih^+i||/II?^., f^H\g^H^+\ 

(A.5) 



Proof. The proof of (2.4) is easily obtained by the calculations of the previous 
proof. To prove (2.5), it is only step (i) of the previous proof which is no longer 
true. However, this is also the reason for the k + 1 on g. That is, step (i) is replaced 
by a simpler Holder inequality. 

□ 
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